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1 Introduction 



The 't Hooft model is a name given to the QCD in two dimensions with 
infinite number of colors [TJ. Among many features of the theory we will 
be interested in the following one: it is shown in pQ that mesons and their 
spectrum can be exactly described in terms of solutions of a certain integral 
equation. Although it can be be solved numerically [21 El lU E], this equa- 
tion, due its elegance and simplicity, deserves an analytical study One is 
interested in analytic properties of the spectrum as functions of parameters 
of theory, that is, the coupling constant and masses of quarks, taken as com- 
plex variables. With this goal in mind, the authors of [H] proposed a novel 
method for studying analytic properties of solutions of the equation, more 
specifically the spectrum of the problem. Although they consider the case 
when quark masses have special values, the method they developed is general 
and many conclusions about properties of the spectrum are applicable or can 
be extended to a general case without any obvious difficulties. 

't Hooft's theory for mesons is not a unique example of theories with 
similar properties. For example, in [BJ [5] a Bethe-Salpeter equation for the 
mesons of the Ising field theory (IFT) was proposed. From the particle 
content point of view, theories are not very different, each being a two- 
dimensional theory of quarks with a confining interaction. At the same time, 
whereas the bound state equation in the IFT is written in a certain, rather 
limiting approximation, the equations that appear in 't Hooft's theory are 
exact. This fact, together with results of studies of 't Hooft's theory, might 
give us some intuition as to what to expect in the full IFT. 

This work can be seen as a preliminary step in studying the case of quarks 
with generic masses. Here we restrict ourselves to the case when quarks have 
equal yet arbitrary masses. The goal of the paper is to explore asymptotic 
properties of the spectrum. For this, we develop low-energy and semiclassical 
expansions for the spectrum, study their properties and compare results with 
numerical analysis. 

This paper is organized as follows. First, we formulate the problem. Next 
section is devoted to studying the weak-coupling limit, which, depending on 
certain conditions, leads either to the low-energy or semiclassical expansions. 
Finally, we compare analytical results with the numerical solution. In the 
appendix we provide another way to derive the low-energy expansion. 
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2 Set-up 



't Hooft's model, or two-dimensional QCD with gauge group SU(N), is de- 
fined by the following Lagrangian 

C = ~ Tr F, U F^ + ^(i-fDn - m { a) )^ a \ (2.1) 

here ijj^ are quarks with bare masses , different in general; the field 
strength F^ u = d^A v — d v A fl + i[A^, A u ] and the covariant derivative = 
+ iA^ are defined in terms of the gauge potential A^ given hy N x N 
hermitian traceless matrices. In the large N limit, 't Hooft derived the 
Bethe-Salpeter equation for mesons, two-quark bound states of the theory [TJ. 
This equation, though originally obtained in the light-cone gauge, is actu- 
ally gauge-invariant since it describes gauge-invariant objects. When quarks 
have equal masses rn^ = m , which is the case we will consider, 't Hooft's 
equation has the form 

MVM = ^4 - 9 -f (2.2) 
x(l-x) ttJq {y-x} 2 

with ip(x) being the wave function of the meson with the mass M and m 2 = 
m o ~ Q 2 / 71 standing for the renormalized quark mass. It is this equation, 
Eq. ( 12. 2p . that is the main object of our analysis. Written in terms of 
dimensionless quantities, a = M 2 /4m 2 and A = g 2 /2m 2 , the equation ( 12.2ft 
becomes 

M*) = " 4 l / {V \^ . (2.3) 

yi ; Ax{l-x) Jo (y-x) 2 2n 1 ; 

The bound-state equation can be viewed as an eigenvalue problem Hip = 
oap for a suitably defined hamiltonian H acting in a certain Hilbert space. 
Solutions of the problem can be classified according to the symmetry x — > 
1 — x of the equation. 

The overall goal of this work is to provide some insight into analytical 
properties of the spectrum, that is, properties of a as a function of A taken 
complex. As we have pointed out in the introduction, at the current stage 
we will be mainly interested in asymptotic properties of the spectrum leaving 
other, more interesting questions for later. 
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3 Weak-coupling expansion 

In the weak-coupling limit A — > 0, each meson's mass M n will approach 2m 
from above, i.e., a approaches l+o- Depending on how A goes to 0, one 
gets different expansions. For example, if n is kept fixed, this limit will give 
us the low-energy expansion. Such expansion is known to be in t = A 1 / 3 . 
At the same time if one looks at high levels, with n > 1/A, one obtains 
the semiclassical expansion giving corrections to the already known Bohr- 
Sommerfeld equation for the spectrum. 

In our analysis for consistency we will use notations adopted in [8]. 

It is convenient to go to the rapidity space, Q = \ log tz~. Then Eq.f l2.3p 
takes the form 




(3.1) 



Define the hamiltonian H as 




(3.2) 



then Eq.f l3.ip becomes equivalent to 



H<p(8) = a<p(9) 



(3.3) 



with <f (6) belonging to the Hilbert space with the metric 




(3.4) 



For ip(x) we use the following ansatz 




(3.5) 



in the odd sector, and 




(3.6) 



in the even sector, where 



S(/3) = atanh/3-/3. 



(3.7) 
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A even (9) = / dPe lS ^ x 



Both ansatz functions can be shown to be normalizable with respect to the 
metric (13. 4p . 

In the further analysis we will be working with the following quantities 
sinhfl = (H - a)A(B) (3.8) 

and 

- coshtf A even (6) = (H- a)ipMjO) (3.9) 

Computation along the lines of the one shown in Appendix B of [8] gives 
the following expressions for these functions 

a sinh ft iX cosh 9 (cosh (3 cosh 9 + 1) 

cosh 2 (cosh p + cosh 9) 2 

(3.11) 

We develop the small-A series for ( 13.1 Oft and ( 13. lip . At first few orders, 
terms in these expansions will have no dependence on 9, directly giving the 
spectra a n as zeroes of A odd / even {9) = A odd / even (9\a). But in general, these 
terms will explicitly depend on 9. In such situations one has to proceed as 
follows. Let {(fin} be the complete orthonormal set of eigenvectors of H: 
Hip n = a n ip n . Then the following quantity, set to zero, 

C n;odd/even (a) = ((f n , (H - a)p2u /ew J ( 3 - 12 ) 

will play the role of the spectral condition with a n being its zeroes. For 
small values of A, ^M/even provides a good approximation for ip n . As A 
grows, corrections can be found by iterating the following expression 

(3 - 13) 

Finding (p n order by order and plugging it in (I3.12p . one obtains corrections 
to the spectral condition. 
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3.1 Low-energy expansion 

Assume that a is close to 1 and write it as a = 1 + zt 2 + Ylk=3 e kt k - in this 
case one sees that the main contribution to the integrals (I3.10P and (13. lip 
comes when ~ t. Rescaling = —ut, 9 = —vt, and observing that 

' S! ° - S^ + S^u), (3.14) 



A 



5(/3) 



where So(u) — ^ — zu, hence giving the Airy function, and S\(u) — x 
So(u) = 0(t), one can expand A(#) in t. Subsequent setting terms in each 
order of t equal to zero gives conditions sufficient for determining a. 

After a straightforward computation one obtains low-energy expansions 
for a in both sectors. In the odd sector the first terms will have the form 

where z is such that Ai(— z) = 0. Whereas for even solutions one obtains 
cw = 1 + th + t 4 [?r + 1 V t 6 + 4 - 7^) +■■ (3- 16) 



^5 5zJ V 175 25 502 , 

with z being the solution of Ai'(— z) = 0. 

As above, one should mention that except for the first few terms of the 
expansion of A in t, expressions at powers of t will depend on v , that is on 9. 
In these cases one has to follow the procedure outlined above. We just note 
that such dependence appears for the first time at t 8 -order in the odd sector 
and £ 6 -order in the even. 

Another observation is that as opposed to the IFT, with odd k seem 
to be absent in low orders both in the odd and even cases (at least up to 
k = 15). 

Looking at the coefficients in (13. 15ft and (13. 16ft . one can see the general 
structure of e^, polynomial in z or in z and 1/z, respectively. At higher levels 
z n behave as z n ~ n 2//3 both for odd and even solutions. Hence, the leading 
terms in e^, terms with the largest power of z, which come from the term at 
A in the brackets in (I3.10P and (13. lip , depend on the combination [n l ^t) 2m . 
Moreover, one can notice that numerical coefficients are the same for both 
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odd and even sectors. This is expected because the sum of these leading 
terms is given by the Bohr-Sommerfeld formula, the leading contribution in 
the semiclassical expansion. Therefore, one would expect that summing sub- 
leading terms should reproduce contributions of the semiclassical expansion 
of higher orders. 



3.2 Semiclassics 

When A goes to zero and a is not necessarily close to 1, one obtains the semi- 
classical approximation. Technically it is reduced to calculating the integrals 
( 13.1 Op and (13. lip using the saddle-point method. 
In both sectors the integral for Aiff) has the form 

/■oo 

A(0) = / d[3f(P\9)e iSW/x . (3.17) 



There are two saddle points = ±i? with •& being the positive solution of 

a = cosh 2 ??. (3.18) 

With A±(8) being contributions from ±i9, the integral can be rewritten as 

A(0) = A + (e) + A46). (3.19) 

In each case, odd or even, one has to compute only A + (6), whereas A_(6) = 
±A* h (9), here +/— is in the odd and even sectors, respectively. 

The leading order calculation yields the following expressions in the odd 
sector 

Ao ^ )v/ ^-cos(^)/A-vr/4), (3.20) 



2VttA cosh 3 i? 
and in the even 

A even (0) 



2iV 7rA sinh d cosh d 

with 



sin (£(tf)/A - tt/4) (3.21) 



m = S{P)\p^ = n • ( 3 - 22 ) 



sinh 2d - 2d 
Y 

At this level, there is no dependence on 6, and the spectral condition in the 
leading order will have the Bohr-Sommerfeld form, S($ n ) = n\(n — no), here 
n = 1, 2, . . ., tiq is equal either to 1/4 for the odd sector or 3/4 for the even. 



6 



In contrast to the IFT, ^-dependence becomes relevant already at the 
one-loop level. Computations at this order give the following expression for 
the odd case 

A " w( ^ vsmh >} : cos (S($)/\ - 7r/4 + \A odd (9\#) + ASijowO?)) , (3-23) 



2v / vrAcosh 2 ^ 
where 



cosh 2 9 tanh -d 
= - (coahg + coBh^ ' (3 ' 24) 



SW*) = + Ao \ 3q (-9 tanh 4 + 6 tanh 2 # - 5) . (3.25) 

4 48 tanh 3 1? v ' 

In the even sector the expression is 

^even [9) 

2i\J 7rA sinh d cosh 9 

= sin (S(tf)/A - tt/4 + AA evm (0|tf) + AS^O?)) (3.26) 

with 

A eve J9m = coshe ( coshdcosh $+ 1 ) ( 3 27 ) 

sinh i? (cosh 9 + cosh i?) 2 ' 

Si-eveJti) = ~^r~ + —5- (-9 tanh 4 •& + 6 tanh 2 - 5) . (3.28) 

' V ; 4tanh 3 tf 48 tanh 3 $ K ' 

In this situation one again has to use the method based on (13.121) . At this 

order, when ip n is approximated by ^ d \/ evenl the integral in (13.121) at small 

A will have the main contribution coming from the saddle points 9 = ±i? 

giving the following quantization conditions in the odd sector 

^ H» = vrA(n-l/4)-A 2 {A odd {d n \d n ) + S 1;O(M (0 n )) + O(A 3 ) (3.29) 

and in the even sector 
sinh 2$ n - 2$ n 

1 n - = n\(n-3/4)- 

- A 2 (A even (# n \K) + SWM) + 0(A 3 ). (3.30) 

Once fin is found, the eigenvalues a n are determined using (I3.18p . 

Calculation in higher orders is technically similar, albeit more involved. 

In principle, using ( 13 .29|) and (I3.30P one can re-derive low-energy expan- 
sions (I3.15P and (I3.16p . respectively. 
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4 Numerical results 

Over the years there have been proposed different approaches for numerical 
solution of the problem, each giving solution to 't Hooft's equation with 
any degree of accuracy. In [8] another method, based on discretizing the 
Fourier-transformed version of 't Hooft's equation, was suggested. Applying 
the Fourier transform to Eq. fl3.ll) 



with K{y) = z//(2sinh^). The main feature of this version of 't Hooft's 
equation is the smoothness of its kernel allowing to discretize the equation 
directly. Even though this method can be not as efficient as others, its ease 
makes up for it. One should mention that this is not the first time this 
parametrization has been used. It has already appeared, for example, in [7] 
and recently in [9]. 

In Tables 1 and 2 we list values for a, i.e., for M 2 in units of 4m 2 , for first 
few levels, both even and odd, at different values of A. Numerical results 
are obtained by discretizing Eq. ( 14.2f) with the number of steps N = 3000 
and step size e = 0.1. This choice of discretization parameters provides 
more than sufficient accuracy in this range of A, at least for the first levels. 
They are compared to the results obtained by means of derived expansions 
(I3.15P and (13.161) . truncated at A 10 / 3 in the odd sector and at A 8//3 in the even 
sector. One can see that within the given accuracy numerical results match 
the expected values with deviations determined by the order of the terms 
dropped from the low-energy expansion. 

It is worthwhile to see how the expansions (I3.15P and (13.161) approach 
the numerical values. For the lowest odd and even levels, which we choose 
because zeroes of Ai(— z) and Ai'(— z) are smallest for them thus making 
specific features of the expansions, especially in the even case, more promi- 
nent, this can be observed in Tables 3 and 4, respectively. Here a[^ et)efi 
denotes estimates obtained using (13. 15[) and (I3.16P by keeping terms up to 
t k included. 




(4.1) 



one obtains 




(4.2) 
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Table 1: a n for first odd levels. The first value is numerical; the second is 
obtained using the truncated low-energy expansion (13.151) ; the third value is 
the order of the difference between the two. 



A 


0.0001 


0.001 


0.01 


0.1 


n = 1 


1.00504237411651 
1.00504237411666 
lO" 13 


1.02349036201 

1.02349036216 
10 -io 


1.11087634 
1.11087649 

io- 7 


1.55404 
1.55418 
IO" 4 


n = 2 


1.00882272353899 
1.00882272353917 
lO" 13 


1.04121272921 

1.04121272942 
10 -io 


1.19684999 
1.19685027 

io- 7 


2.02721 
2.02759 
IO" 4 


n = 3 


1.01192195048315 
1.01192195048344 
lO" 13 


1.05581244361 

1.05581244407 
10 -io 


1.26911411 
1.26911489 
10~ 6 


2.44971 
2.45097 
IO" 4 


n = 4 


1.01466422536562 
1.01466422536618 
lO" 12 


1.06878314251 
1.06878314354 

io- 9 


1.33436550 
1.33436743 
IO" 6 


2.84797 
2.85117 
10~ 3 


n = 5 


1.01717361844767 
1.01717361844723 

io- 13 


1.08069521197 
1.08069521405 

io- 9 


1.39513452 
1.39513856 
IO" 6 


3.23152 
3.23815 

IO" 2 



Tables 5 and 6 contain values of a n obtained through semiclassical expan- 
sions, (I3.29P and (I3.30p . As expected, with n increasing, semiclassical results 
approach values obtained through discretization. 

5 Conclusion 

In this work we apply methods proposed in [61 18| to the study of 't Hooft's 
model. In particular, we study asymptotic properties of the spectrum of 
mesons and obtain the low-energy and semiclassical expansions for meson 
masses. We compare these results to the numerical solution of the problem. 

This study has to be considered as a part of an extended project. Methods 
used in this work, together with the approach proposed in [9], will likely 
provide a way to tackle the general case of the problem with arbitrary quark 
masses. The project itself has many possible ramifications. Some of them 
are related to questions about analytical properties of the spectra which are 
preserved if one considers the case of number of colors large but finite. In such 
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Table 2: Same as in Table 1, but for first even levels. 



A 


0.0001 


0.001 


0.01 


0.1 


n = 1 


1.002196798442 

1 nmi QR7QSJ/1 81 

i.uuzryo i yo4oi 

io- n 


1.0102283943 
1 nimosj/inop; 

1.U1UZZO4UZ0 

IO" 8 


1.0481649 

l.U4olDDD 

IO" 6 


1.2387 
i.zoyr 
10" 4 




1.007008105144 
i nn7nnsi dm /LQ 

i.uu i uuoiuji^y 
lO" 12 


1.0326986818 

i.uoiuyouou i 
1Q -9 


1.1554011 

IO" 6 


1.7964 

IO" 3 


n = 3 


1.010406331508 
1.010406331481 

io- n 


1.0486680231 
1.0486680086 

io- 8 


1.2336564 
1.2336495 
IO" 5 


2.2408 
2.2381 
10" 3 


n = 4 


1.013313818353 

1.013313818267 
10 -io 


1.0623921932 
1.0623921521 

io- 8 


1.3021422 
1.3021238 
10~ 5 


2.6502 
2.6430 
10- 2 


n = 5 


1.015933386177 

1.015933385994 
10 -io 


1.0748048015 
1.0748047167 

io- 7 


1.3650255 
1.3649880 
IO" 5 


3.0406 
3.0266 
IO -2 



cases one is primarily interested in the nature of singularities. Others deal 
with possible generalizations of these methods to more general field theories. 
For example, one can consider theories with vacuum degeneracy, possessing 
certain integrabilty which becomes broken by perturbations that typically 
give rise to the confining force between particles of original theories. IFT 
and 't Hooft's theory provide sufficiently simple examples of such theories. 
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A Coordinate-space representation 

Dependence of the terms in the low-energy expansion on roots of Ai'(— z) in 
( 13.161) looks drastically different from the behavior in the odd sector. There- 
fore, it is worthwhile to provide another way to reproduce it. 
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Table 3: Step-by-step approximation of the meson mass for the the lowest 
odd level obtained by successively including higher terms in (13.1 5ft ; the last 
value is the numerical result. 



A 


0.0001 


0.001 


0.01 


0.1 


(2) 
a odd 


1.00503729971412 


1.02338107410 


1.10852533 


1.50373 


(4) 
a odd 


1.00504237459180 


1.02349040903 


1.11088088 


1.55448 


(6) 
a odd 


1.00504237411491 


1.02349036134 


1.11087611 


1.55400 


a {8) 

"odd 


1.00504237411666 


1.02349036215 


1.11087649 


1.55418 


(num) 
a odd 


1.00504237411651 


1.02349036201 


1.11087634 


1.55404 



Table 4: Same as in Table 3, but for the lowest even level. 



A 


0.0001 


0.001 


0.01 


0.1 


(2) 


1.002194922920 


1.010187930 


1.0472882 


1.21949 


(4) 
eleven 


1.002196797651 


1.010228320 


1.0481584 


1.23824 


a (6) 

"even 


1.002196798481 


1.010228403 


1.0481666 


1.23907 


(num) 
eleven 


1.002196798442 


1.010228394 


1.0481649 


1.23870 



In [6] a coordinate-space formulation for the bound-state equation was 
used to obtain the low-energy series for the spectrum of the IFT. Similarly 
we map 't Hooft's equation into the coordinate space. For simplicity we 
choose Eq. fl2.3p as a starting point for our analysis. In u = 2x — 1, this 
equation becomes 

a<p(u) = f±-^-2\f ^ ' — . A.l 

1 — U z J_ x [U — V) 2 2,71 

Rescaling u — tk, v — tp, setting A = t 3 , expanding everything in t, and 
rearranging the terms, one gets 

(1 _ a + + t A k A + ... ) V (tk) = 2t 2 T T^%^. (A.2) 



(p- £;) 2 27r' 

Here we changed the limits of integration from (— 1/t, 1/t) to (— oo, oo), an 
operation legitimate at 0(t 4 ). 

Now this equation in the configuration space becomes 
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Table 5: a„ for first odd levels from semiclassics. 



A 


0.0001 


0.001 


0.01 


0.1 


n = 


1 


1.005045002 


1.023502645 


1.110935189 


1.554453706 


n = 


2 


1.008822931 


1.041213704 


1.196854870 


2.027350588 


n = 


3 


1.011921999 


1.055812673 


1.269115381 


2.449829823 


n = 


4 


1.014664243 


1.068783226 


1.334366056 


2.848080751 


n = 


5 


1.017173627 


1.080695251 


1.395134852 


3.231630623 



Table 6: a„ for first even levels from semiclassics. 



A 


0.0001 


0.001 


0.01 


0.1 


n = 


1 


1.002079962 


1.009690485 


1.045761963 


1.229713753 


n = 


2 


1.007007280 


1.032694849 


1.155383793 


1.796538843 


n = 


3 


1.010406206 


1.048667441 


1.233654091 


2.240958482 


n = 


4 


1.013313781 


1.062392019 


1.302141769 


2.650298641 


n = 


5 


1.015933371 


1.074804730 


1.365025460 


3.040735165 



Here the Fourier transform conventions in terms of the rescaled distance 
between quarks x = tx rea [ are 

dk 



ip(x) 



2vr 



e lkx <p(tk), 



and we let 1 — a = —et 2 . 

This equation flA.31) has a symmetry x — > —x, hence its spectrum can 
divided into even and odd parts. 

To find a perturbative solution for it, one first solves the auxiliary equa- 
tion 

d 2 o d 4 A d % 







(A.4) 



with 



F(x) = F (x) + t 2 Fi(x) + t i F 2 {x) + . 
The bounded perturbative solution for (IA.4|) is 

t 2 

F(x) = Ai(x) H UxAi(x) +x 2 Ai'(x)) + 

5 



4 f fx 5 5x 2 \ w , /9a; 3 6 \ A ... , 

+ 1 4 h — M(x) + — + — Ai'fx 

VV50 7 J V ' V 35 35/ V ' 



(A.5) 
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The following two claims provide both the solutions for (1A.3j) and the 
spectra for even and odd sectors. In each case we assume that F(x) is the 
perturbative solution (IA.5j) for (IA.4j) . 



Claim 1. ipodd(x) = sgax F(\x\ — e) provides the odd perturbative 
solution of flA.3j) if F(x) is such that 



F(-e) = 0(t 2 ), 

F{-e) - t 2 F"(-e) = 0(t 4 ) } (A.6) 
F(-e) - t 2 F"{-e) + t 4 F w (-e) = 0(t 6 ). 



Claim 2. ip e ven{.%) = F(\x\ —e) is the even perturbative solution of (1A.3j) 
if F(x) satisfies 

F'(-e)=0(t 2 ), 

F'(-e) - t 2 F"'(-e) = 0(t 4 ), (A.7) 
F\-e) - t 2 F"'{-e) + t 4 F v (-e) = 0(t 6 ). 

The proof of each claim is straightforward. Higher order conditions can 
be easily written but they are not necessary at the order we are interested 
in. 

In both cases conditions on F, (IA.6I) and flA.7j) . are sufficient for deter- 
mining the spectra. It is easy to check that using them one obtains the 
following expansion of a in the odd sector 

>2„ i ^ 2 x Jo ( ^ 3 6 



a odd = 1 + t 2 z + — + t b ( -y^z + 35 J + • • ■ (A ' 8) 



and in the even sector 



C'l . /- I / ^ 2 i ^ i j-6 / ^ 3 3 



where z is the solution of Ai(— z) = and Ai'(-z) = 0, respectively, thus 
matching low-energy expansions (13.151) and (13. 16[) derived above at the chosen 
order in t. 
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